Velocity Fluctuations in Electrostatically Driven Granular Media 
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We study experimentally the particle velocity fluctuations in an electrostatically driven dilute 
granular gas. The velocity distribution functions have strong deviations from Maxwellian form in a 
wide range of parameters. We have found that the tails of the distribution functions are consistent 
with a stretched exponential law with typical exponents of the order 3/2. Molecular dynamic 
simulations shows qualitative agreement with experimental data. Our results suggest that this non- 
Gaussian behavior is typical for most inelastic gases with both short and long range interactions. 
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Despite extensive study over the preceding decade, a 
fundamental understanding of the dynamics of granu- 
lar materials still poses a challenge for physicists and 
engineers Driven granular materials exhibit com- 

plex behavior that resemble some aspects of conventional 
solids, liquids, and gases, yet also demonstrate some con- 
siderable differences 

Recent experimental studies with vibrationally-driven 
granular gases revealed surprising deviations in the par- 
ticle distribution function from the Maxwell distribution 
law (5-10 . These deviations were attributed to the ef- 
fects of dissipation due to inelasticity of inter-particle col- 
lisions. In particular, Ref. 0] reported an "universal" 
stretched exponential law 



P(v) ~ exp[-\v/v Q \ ( 



(1) 



with the exponent £ = 3/2. Here v is particle veloc- 
ity, P is the velocity distribution function and vq is the 
"thermal velocity" . This behavior was observed in a 
wide range of frequencies and amplitudes of vibration 
and for different densities of the granular gas and agrees 
with the theoretical prediction of Ref. JTlj ] . Deviations 
from Maxwellian behavior also has been observed exper- 
imentally in different geometries and for different driving 
conditions in Refs. and in numerical simulations 

fl2| llif . These studies suggest that the deviations from 
a Maxwell distribution are the result of short-range in- 
elastic hard-core collisions between particles. 

Interactions between particles often are not reduced 
to simple hard-core collisions. Fascinating collective be- 
havior appears when small particles acquire an electric 
charge and respond to competing long-range electromag- 
netic and short range contact forces. An important ques- 
tion is whether non-Maxwellian distributions are typical 
only for granular gases with hard-core collisions or could 
it also be observed for general dissipative gases with both 
short-range and long-range interactions. 

The electrostatic excitation of granular media offers 
unique new opportunities compared to traditional vibra- 
tion techniques which have been developed to explore 
granular dynamics. It enables one to deal with extremely 



fine powders which are not easily controlled by mechani- 
cal methods. Electrostatic driving makes use of these bulk 
forces, and allows control of the ratio between long-range 
electric forces and short-range collisions by changing the 
amplitude and the frequency of the applied electric field. 
Our previous studies with electrostatically-driven granu- 
lar media revealed a hysteretic phase transition from the 
immobile condensed state (granular solid) to a fluidized 
dilated state (granular gas) with a changing applied elec- 
tric field |l5|] . A spontaneous precipitation of dense clus- 
ters from the gas phase and a subsequent coarsening - 
coagulation of these clusters - is observed in a certain 
region of the electric field values. The strong effect of 
humidity on dynamics of electrostatically-driven granu- 
lar materials was studied in Ref. | p| . 

In this Letter we study the particle velocity distribu- 
tions in electrostatically driven granular media. We have 
found that in a wide range of parameters the particle 
velocity distribution function is strongly non-Maxwellian 
and is well-approximated by the stretched exponential 
law P(v) ~ exp(— \v/vq\ 3 / 2 )- We performed molecu- 
lar dynamics simulations of conducting particles in an 
ac electric field and have obtained qualitative agreement 
with the experiment. We conclude from our results that 
the tails of the velocity distributions for driven granular 
gases in general exhibit non-Maxwellian behavior and are 
not limited to the systems with just hard-core collisions. 

Our experimental setup is similar to that in Ref. 
iPlpll . Particles are placed between the plates of a 
large capacitor which is energized by a constant (dc) or 
alternating (ac) electric field E = Eocos(u>t), see Fig. 
[|. To provide optical access to the cell, the capacitor 
plates were made of glass with a clear conductive coat- 
ing. We used 11 x 11 cm capacitor plates with a spacing 
of 1.5 mm. The particles consisted of 165 y,m conducting 
bronze spheres. The field amplitude Eq varied from to 
10 kV/cm and the frequencies / = ui/2ir on the interval 
of to 120 Hz. The total number of particles in the cell 
varied between 10 5 and 10 6 . 
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momentum transfer occurs often without actual collision. 
While there are hard sphere collisions in the system as 
well, this interaction occurs at a center to center distance 
of approximately two particle diameters. 



Transparent Electrode 



E=E„cos(2Kft) 



FIG. 1. Block-diagram of experimental apparatus. 

Conducting particles acquire a surface charge when 
they are in contact with the capacitor plate. As the 
magnitude of the electric field in the capacitor exceeds 
a critical value, Ex, the upward electric force overcomes 
gravity mg (m is the mass of the particle, g is the grav- 
ity) and lifts the charged particles. When grains hit the 
upper plate, they deposit their charge and fall back. By 
applying an alternating electric field E = Eq sin(27rf), 
and adjusting its frequency /, one can control the ver- 
tical excursion of particles by effectively turning them 
back before they collide with the upper plate. Thus, by 
increasing the frequency of the electric field, the parti- 
cles can be confined in a relatively thin layer, i.e. the 
granular gas is quasi-two-dimensional at high frequencies 
and three-dimensional at low frequencies. It effectively 
allows for control over the number of collisions and, there- 
fore, the contributions from long-range and short-range 
interactions. This control over the vertical extent of the 
motion of the particles also changes the manner in which 
energy is transferred from the vertical to the horizontal 
through collisions. From the comparison of typical ki- 
netic energy acquired by the particle and electric energy 
of inter-particle interaction it is possible to show that the 
long-range electric interaction becomes dominant with 
respect to short-range collisions at the frequencies larger 
than fo « 100 Hz for given particle size. 

In our experiment we extracted horizontal particles ve- 
locities via high-speed image analysis. Pictures were ob- 
tained at the rate 1000 frames per second from a camera 
mounted to a microscope suspended vertically above the 
cell and particles positions were resolved to sub-pixel res- 
olution. Inter-particle and particle-boundary collisions 
that introduce sudden changes in momenta of the par- 
ticles were filtered from the distributions in a manner 
similar to Ref ||@]. These events were only a small frac- 
tion of the total number of measurements due to the low 
density at which they were acquired. An ensemble aver- 
age for each of the velocity distributions were obtained 
with nearly 10 6 data points. 

The effect of long-range electric forces is illustrated in 
Fig. 0. In contrast to vibrationally-driven systems, the 




FIG. 2. Composite images of two particles collision created 
by overlaying several sequential images taken with the high 
speed camera for a driving frequency of 50Hz. Each of the 
three images demonstrates a different collision between two 
particles that occurs at a distance larger than a ball diam- 
eter. At this relatively high frequency the particles have a 
small vertical excursion and the focal plane of the microscope 
is set between the two plates of the cell. The particles are 
completely illuminated which results in the smudged streaks 
of varying intensity in the particle tracks. 

A summary of experimental results is presented in 
Figs. U and |]. We find for various values of frequency 
/ and electric field amplitude Eq the normalized veloc- 
ity distributions both in logarithmic and linear scales lay 
practically on a top of each other, see Fig. 0. There is 
a very small difference in granular temperature T g be- 
tween x and y directions, which is probably related to 
some small tilt of the cell with respect to gravity. There 
is an increase in the flatness of the velocity distributions, 
F = (v 4 ) / (v 2 ) 2 — 3, as the driving frequency is decreased 
at a fixed driving amplitude, see Fig. ||,inset. (As defined 
here, a flatness of corresponds to a Gaussian distribu- 
tion). In a mechanically shaken granular layer, an in- 
crease in flatness was observed as the shaking amplitude 
was decreased at constant frequency, which corresponded 
to a lower horizontal granular temperature Q). In the 
electrostatically driven system, the increase in flatness 
occurs as the horizontal granular temperature increases 
as there is more kinetic energy to transfer from the hori- 
zontal to the vertical direction because of the increase in 
vertical excursion at lower frequencies. 



2 



-1 



0<3 



o 



0.0 



• 


0.6 






if 




0.4 


* # 










0.2 


# # 





v/v n 



<> 



4- 



tion was built from 2,500,000 sample points. 
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FIG. 3. Main plot: The logarithms of horizontal veloc- 
ity distribution P(v) vs normalized velocity v/vo, where 



correspond to / = 50 Hz and Eo = 8.83 kV/cm, stars to 
/ = 120 Hz and E = 8.83 kV/cm, diamonds to / = 45 Hz 
and Eo = 7.3 kV/cm. Dashed line show the best fit accord- 
ing to Eq. (Q) with £ ~ 1.51. Inset: The horizontal velocity 
distribution P(v) vs normalized velocity v/vo in linear scale. 

Fig. |] shows the density distribution function g(r) = 
(p(r)/p(0))/(p) 2 for different values of frequency. The 
volume density for the experimental conditions is ex- 
tremely dilute (w 7%) as can be seen by the lack of 
particle correlations with the exception of the excluded 
volume of the finite particles. Consequently, the func- 
tion g(r) has almost step-like form for low frequencies. 
One can also see that with the increase of frequency the 
contribution from long-range interactions also increases. 
This fact is manifested in the smearing of the step-like 
structure of g(r), which can be interpreted as an increase 
of the effective crossection for particle interactions. 

In order to obtain an insight into the problem we 
performed molecular dynamics simulations of conduct- 
ing particles in an applied ac electric field. We modified 
the MD code used in our earlier paper fl5]| . In particular, 
we introduced a finite roughness of the bottom plate and 
finite restitution coefficient r between the particles and 
particles and walls (r = 0.8 and r — 0.6 correspondingly) 
ipTf . We changed the restitution coefficient and rough- 
ness over a wide range, but the results appear to be very 
robust to the specific choice of r. We simulated up to 
1000 particles in the cell 80 x 80 particles diameters in 
x — y plane, electric field amplitude Eq was 2.6i?i, where 
Ei is the first critical field value (for E = E\ upward 
eclectic force acting on isolated particle overcomes grav- 
ity, see The velocity is scaled by v = ^fga and time 
by t — \/ a/g, where a is the radius of particle. Typical 
simulation time was 10,000 i, and the velocity distribu- 
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FIG. 4. The density distribution function g(r) for different 
values of /. Inset: "thermal velocity" vq (solid line) and 
flatness F (dashed line) vs / for Eo = 8.83 kV/cm. 

Selected results are presented in Fig. ||. In the wide 
range of parameters u>, Eq we have found non-Maxwellian 
velocity distributions. For relatively high frequencies 
lj = 0.5 (corresponding approximately to / = 30 Hz for 
our experimental parameters) the velocity distribution is 
indeed approximated by Eq. (Q) with ( close to 3/2, as 
in the experiment. For a smaller frequency (lo = 0.2) 
we have found that P(v) is consistent with a simple ex- 
ponential law: P(v) ~ exp(— v/vq), i.e. £ = 1. It is 
plausible that the crossover to exponent £ = 3/2 occurs 
for the high velocity tail and would explain the gentle in- 
crease in flatness as the frequency is reduced. However, 
one likely needs substantially longer simulations in order 
to resolve this crossover. 

Our simulation results for "high frequency" yield an 
exponent C ~ 1-5 and are in a good agreement with 
experiment. Unfortunately, we cannot verify the "low 
frequency" exponent ( = 1 because our current exper- 
imental apparatus does not allow us to extract reliable 
particles positions for the frequencies / below 30 Hz. For 
small frequencies the vertical particle displacements are 
much higher than the particle size. As a result, the par- 
ticles go out of the focal plane of the microscope and it 
is impossible to sufficiently resolve horizontal positions. 

Let us briefly discuss the results. Universal high ve- 
locity tails of the distribution Pcxp(— \v/vq\ 3 ^ 2 ) were 
derived by van Noije and Ernst |ll| from the Enskog- 
Boltzmann equation for the uniformly heated state of a 
2D weakly-inelastic granular gas. According to pH , in 
the steady state, the averaged distribution function sat- 
isfies the equation 
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Dd 2 v P = -I(P,P) (2) 

where I{P, P) is the pair collision integral and D is pro- 
portional to the strength of external "thermal noise" . In 
the large velocity limit Ref. derived I(P, P) ~ vP 
since the fast particle typically collides with particles in 
the "thermal range" , which immediately results in Eq. 
(|l|) for the high velocity tail. 

Although there is no "thermal noise" in our exper- 
iment, multiple collisions between particles, the finite 
roughness of the plates, and long-range interactions re- 
sult in sufficient randomization of particle velocities and 
presumably have an overall effect similar to thermal 
noise. The argument that the collision integral for high 
velocities is reduced to vP appears to be very robust and 
does not depend on the specifics of the model. Thus, we 
can expect the above stretched exponential velocity tails 
to be characteristic of our system as well. 



is well-approximated by the stretched exponential law 
P{v) ~ exp(— |ii/i;o| 3 / 2 ). Our results suggest that such 
velocity distributions should be generic for dissipative 
gases with both short and long range interactions and 
may possibly have some relevance for completely differ- 
ent systems, such as dusty plasmas, colloids, etc. 
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FIG. 5. The logarithms of horizontal velocity distribu- 
tion P(v x ) (i-component only) for two different values of fre- 
quency uj. Solid lines show numerical results, dashed lines 
are fits \ogP(v) — ao — div , with £ « 1.43 for oj — 0.5 and 
C « 1.08 for lj = 0.2. 

The distribution P ~ exp(— v/vq) was derived for a 
freely evolving granular gas |Q . One can speculate that 
our results for the small frequency have some fingerprints 
of a "freely evolving gas" because the particles spend sig- 
nificant time in the flights between plates and exhibit 
some sort of "free cooling" , whereas for the higher fre- 
quencies the particles are confined in a thin layer and 
collisions with the bottom plate occur more often. 

In conclusion, we experimentally and numerically stud- 
ied particle velocity distributions in an electrostatically 
driven granular media. We have found that in a wide 
range of experimental parameters the particle veloc- 
ity distribution function is strongly non-Maxwellian and 
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